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Abstract: Partial sequences of the D-loop and the complete sequences of cytochrome b gene (1 140 bp) of the slow 
lorises (genus Nycticebus ) were undertaken to investigate evolutionary relationships among species of Nycticebus. Se- 
quence analysis results consistently provide new taxonomy evidence at the DNA level for supporting Ratajszczak and 
Groves’ viewpoint that N . intermedus is merely the adult of N . pygmaeus (Ratajszczak, 1998; Groves, 1971). Phyloge- 
netic analysis was performed by means of the combined data and these two separate sequences data, respectively, by us- 
ing various methods, supporting the same topology, in which genus Nycticebus was formed of two clusters. The first clus- 
ter was composed of N . pygmaeus , and the second cluster of N . coucang. It also could provide a new molecular genetic 


evidence to support the view that the genus comprises two species: N.coucang and N . pygmaeus . 
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The slow lorises (genus Nycticebus ) are small 
nocturnal prosimians of family Lorisidae. They belong 
to extant suborder of Strepsirrhine, Primates. The 
slow loris ( N. coucang ) is widely distributed across 
the mainland and islands of Southeast Asia ( Rasmussen 
& Nekaris, 1998). The lesser or pygmy slow loris 
(N . pygmaeus ) is mainly found in Vietnam, Laos, 
China, and possibly also in Cambodia ( Nowak, 
1991, Dang, 1998). The number of the slow lorises 
has severely decreased due to habitat loss and hunting 
(Jurke et al, 1998; Su et al, 1998). N. pygmaeus 
is listed as vulnerable by IUCN (The World Conserva- 
tion Union), itis also listed as one of the first protect- 
ed animals in China. 

There is still controversy over taxonomy in genus 
Nycticebus. Some morphologists support the view that 
there are one to three species in genus Nycticebus 
( Dao, 1960; Groves, 1971; Petter & Petter- 
Rousseaux, 1979; Hill, 1953; Nowak 1991; Su et 
al, 1998). Based on morphological information, the 
slow loris has been described by these as one species, 
N . coucang ( Weisenseel et al, 1998; Hil, 1953; 
Petter & Petter-Rousseaux, 1979). Bonhote de- 
scribed N .pygmaeus as a separate species in 1907 
(Ratajszczak, 1998). Dao (1960) described another 
Indochineses loris species, N .intermedius, which 
was supposed to be intermediate between the very small 
form, N.pygmaeus and the large form, N .coucang 
( Ratajszczak, 1998), although, Groves ( 1971 ) 
claimed that it is merely the adult of N . pygmaeus. 
Groves’ view suggested that Nycticebus comprises two 
species, N.coucang and N.pygmaeus ( Groves, 
1971). Genetic studies from karyotyping, chromoso- 
mal polymorphism data (Warter & Rumpler, 1998; 
Chen et al, 1993), from protein polymorphism data 
(Su et al, 1998), from mitochondrial DNA-RELP 


(Zhang et al, 1993) and restriction maps of ribosomal 
DNA repetypes (Wang et al, 1996), provide evi- 
dence for supporting the view that the genus comprises 
two species, N.coucang and N . pygmaeus . 

Mitochondrial DNA (mtDNA) can be used as a 
good phylogenetic marker in the context of mammalian 
evolution (Saccone et al, 2000) because animal mtD- 
NA is a small molecule that lacks recombination and is 
generally materially inherited, and the rate of animal 
mtDNA evolution is faster than that of single copy nu- 
clear DNA (Fowler et al, 2000; Brown et al, 
1979). Lu et al (2001) reported cytochrome b gene 
sequences from 3 individuals of the slow loris. In this 
study, the partial sequences of D-loop and complete 
sequences of cytochrome b gene of mtDNA were se- 
quenced to investigate evolutionary relationships among 
species of the slow lorises. The molecular data not on- 
ly provides new taxonomic evidence, but should also 
be useful to conservation and genetic management of 
the endangered animal (Su et al, 2001). 


1 Materials and Methods 


1.1 Samples and DNA extraction 

Examined samples were composed of hair biop- 
sies. Hair specimens of seven individuals in slow loris- 
es were used, of which four were N. pygmaeus, one 
was N . intermedius identified according to morphologi- 
cal information, and two were N.coucang. Table 1 
presents species and location of collection of the sam- 
ples. Using Insta Gene Matrix, hair DNA purification 
Kit (Bio-rad Inc.), DNA extraction was performed on 
the hair specimens. The supernatant were kept at — 
20 © for PCR amplification . 
1.2 Sequences from GenBank 

The sequences from GenBank that were used in 
the present study are shown in Table 2. Since all of 


Table 1 Species and localities of collection of the samples in the present study 





Species characterized by morphological information Abbreviated code Locality Origin 
Nycticebus pygmaeus Npyg! KFBG unknown 
Npyg2 KFBG unknown 
Npyg3 KFBG unknown 
Npyg4 KFBG unknown 
N. intermedius Npyg5 KIZ Southeast of Yunnan or Vietnam 
N. coucang Ncou3 UZG Malaysia 
Ncou4 UZG Malaysia 


KFBG: Kadoorie Farm & Botanic Garden Corporation; KIZ: Laboratory of Molecular Evolution and Genome Diversity, Kunming Institute of 
Zoology, the Chinese Academy of Sciences; UZG: Ueno Park, Taito-ku, Tokyo, Japan, Ueno Zoological Gardens. 
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Table 2 Sequences that derived from GenBank were used in the present study 


Species Accession number Abbreviated code 
Nycticebus coucang NC-002765 ncoul 
N. coucang AJ309867 ncou2 
Lemur catta NC-004025 Lemca 
Tarsius bancanus NC-002811 Tarsiban 


them are sequences of mitochondrial genome, both of 
the homology sequences corresponding to the partial D- 
loop sequences and the complete sequences of cy- 
tochrome 6 in the present study were collected from the 
mitochondrial genome. 

1.3 Amplification and sequencing 

Primers L15996 (Yuan et al, 2001) and H16498 
(Ward et al, 1991, Yuan et al, 2001) were em- 
ployed for amplifying the 5’ portion of mtDNA D-loop. 
The complete sequences of cytochrome b were obtained 
by using the universal primers, including two primer 
pairs: L14724 and H15149, and L15136 and H15915 
(Irwin et al, 1991). In order to obtain accurate re- 
sults of cytochrome b, the two fragments amplified by 
those two primer pairs were repetitively sequenced at 
least two times. It was also ensured that those same 
fragments, used to sequence repetitively, were derived 
from the PCR products that were amplified at different 
times. PCR reactions were performed using standard 
buffer conditions, 1.5 mmol/L MgCl,. Reactions were 
subjected to an initial denaturation step at 95 © for 2 
min, followed by 35 cycles at 94 © for 1 min (denatu- 
ration), 42- 50 © for 1 min (annealing) and 72 C 
for | min (extension), and 10 min at 72 © for final 
extension. The different annealing temperatures were 
used with the different samples of slow lorises. 

The PCR products were isolated on a 2%, low 
melting point agorose gel electrophoresis, and purified 
using PCR Preps DNA Purification Kit ( Huashun 
Biolnc. Shanghai). Using BigDye TM Terminator Cy- 
cle Sequencing Kit (Perkin-Elmer), sequencing reac- 
tion was performed (according to the method recom- 
mended by Inc.). Reaction products were purified us- 
ing Sephadex G-50 (ICN Biomedicals Inc. ); resulting 
fragments were analyzed in an ABI Model 377 XL auto- 
mated sequencer (PE Biosystems, Perkin-Elmer, Fos- 
ter City, CA, USA). Each amplified fragment was 
sequenced in both directions. 

1.4 Data analyses 

The partial sequences of mtDNA D-loop and the 
complete sequences of cytochrome b were checked by 
eyes and aligned using DNASTAR package (DNASTAR 


Inc. ). Initial sequence comparisons, base composition 


and measures of variability were performed using MEGA 
version 2.1 (Kumar et al, 2001). The polymorphic 
sites and synonymous/non-synonymous substitutions 
were calculated using software DNASP3.0 (Rozas & 
Rozas, 1999). Phylogenetic analyses were performed 
using PAUP* Version 4.08b ( Swofford, 2001 ). 
Neighbor-joining trees based on P-distance (Nei & Ku- 
mar, 2000) were computed, and 1 000 bootstrap 
replicates were conducted to assess the statistical confi- 
dence of each node. The maximum parsimony analysis 
was performed using heuristic search, and 
the majority - rule consensus tree was obtained from 
1 000 bootstrap replicates. The maximum likelihood 
analysis was also carried out using PAUP. Lemur catta 
(abbreviated as Lemca) and Tarsius bancanus (abbre- 


viated as Tarsiban) were used as outgroups. 
2 Results and Discussion 


2.1 D-Loop sequences analysis 

PCR products of the 5’ portion of mtDNA D-loop, 
adjacent to tRNA-Pro gene, from seven individuals of 
slow lorises, were sequenced and obtained approxi- 
mately 400 bp with an aligned total length. These sev- 
en sequences, together with two corresponding se- 
quences of N . coucang from Genbank, were given se- 
quence analysis. Among the sequences, there were 57 
variable polymorphic sites, and 56 parsimony informa- 
tive polymorphic sites. Only one site was a singleton 
polymorphic site. There was an individual (namely , 
Npyg5), which was identified as N . intermedius, ac- 
cording to morphological information. However, the 
analysis of the D-loop sequences revealed that the se- 
quence of Npyg5 was identical to that of three of four 
individuals of N . pygmaeus (namely, Npygl, Npyg2, 
and Npyg4). They were identified as being the same 
haplotype. This provides new taxonomic evidence with 
molecular genetic study for supporting N . intermedius is 
merely the adult of N . pygmaeus (Ratajszczak, 1998; 
Groves, 1971). 
2.2 Cytochrome b sequences analysis 

PCR products of cytochrome b gene from the same 
seven individuals of slow lorises were sequenced. The 
cytochrome b data set consisted of 1 140 nucleotides for 
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which the codon positions were determined. These se- 
quences were aligned with two sequences of N . coucang 
from GenBank. No insertions or deletions were ob- 
served. No terminal codon was found between the first 
and last codons, and all of these sequences started with 
initial codon * TGA?’ and ended in terminal codon 
‘AGA’. These 1 140 bp of cytochrome b sequences 
coded 379 amino acids in length. It was clear that all 
sequences obtained were for the genuine cytochrome 6 
gene. The analyses of base composition suggested that 
the characteristics are consistent with base composition- 
al patterns found in various mammalian species (Irwin 
et al, 1991). 

Additionally, the analysis of cytochrome b se- 
quences also revealed that the sequence of Npyg5, 
which is designated as N . intermedius according to mor- 
phological information, was identical to that of two in- 
dividuals of NN. pygmaeus Npygl 
Npyg3). They were identified as the same haplotype. 
These results derived from cytochrome 6 sequences 


( namely, and 


were comparable to that of D-loop sequences. 
2.3 Phylogenetic analysis 

GeBank accession numbers of sequences from the 
partial D-loop and the complete cytochrome b gene of 
DNA are AY687887, AY687888, 


mitochondrial 
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Fig. 1 


AY687895 — AY687899 and AY687889 — AY687894, 
AY687900 in seven slow lorises ( Nycticebus ), respec- 
tively. 

The phylogenetic analysis combined partial D-loop 
sequences data with complete sequences of cytochrome 
b gene data was carried out. The PTH test was per- 
formed first. P = 1.0000 was obtained from this test. 
The result revealed that combined data analysis was ap- 
propriate. In the parsimony analysis, there were 311 
parsimony-informative sites found among the nine se- 
quences of genus Nycticebus. The most parsimonious 
tree is shown in Fig.1. It shows no trees that were one 
or two steps less parsimonious. The topology suggested 
that genus Nycticebus was composed of two clusters. 
This was supported by the very high bootstrap value, 
up to 100%. Cluster ĮI was grouped by all individuals 
of N.pygmaeus and one individual of N . intermedius 
(i.e. Npyg5) characterized by morphological informa- 
tion. The bootstrap value was also very high, at 
100% . It is obvious that cluster [ is a monophyletic 
group. Cluster [[ was composed of four sequences of 
N . coucang . This was well supported by high bootstrap 
value, at 100%. Cluster [[ also is a monophyletic 
group. 

The maximum likelihood analysis was also per- 
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The most parsimonious tree of combined partial D-loop sequences data with complete 


cytochrome b gene sequences data of mitochondrial DNA, constructed by the maximum 
parsimony analysis. Numbers on internal branches are bootstrap values derived from 

1 000 replications. Npyg and Ncou are sample code, each number behind them is 

the individual number. ncoul and ncou2 were represented the sequences derived from 
GenBank. Lemur catta (abbreviated as Lemca) and Tarsius bancanus (abbreviated 


as Tarsiban) were used as outgroups. 
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formed. The topology of the maximum likelihood tree 
was consistent with the one of the most parsimonious 
tree. The topology of phylogenetic tree constructed by 
the neighbor-joining method was also accordant with 
that of the most parsimonious tree. 

Whether analysis combining D-loop sequences da- 
ta with complete sequences of cytochrome b gene data 
was performed or analyses were employed each of them 
separately, all of results derived revealed that Nyctice- 
bus is composed of two monophyletic clades that were 
well supported by high bootstrap value. The first clus- 
ter comprised all individuals of N . pygmaeus and one 
individual of N . intermedius; the second cluster was 
grouped with all individuals of N . coucang . The phy- 
logenetic relationship of genus Nycticebus provides new 
molecular genetic evidence with classification of Nyc- 
ticebus. It supports the existence of two different 
species of genus Nycticebus. Groves (1971) suggested 
that Nycticebus consists of two species, N .coucang 
and N. pygmaeus, and this is supported by Ravosa’s 
study (Ravosa, 1998). Genetic studies from kary- 
otyping, chromosomal polymorphism data (Warter & 
Rumpler, 1998; Chen et al, 1993), from protein 
polymorphism data (Su et al, 1998), from mitochon- 
drial DNA-RELP (Zhang et al, 1993) and restriction 
maps of ribosomal DNA repetypes ( Wang et al, 
1996), provide evidences for supporting the view that 
the genus comprises two species: N.coucang and 
N . pygmaeus . Weisenseel et al (1998) first performed 
comparison of all reproductive parameters between 
N. pygmaeus and N.coucang. The results of this 
work support the existence of two different species of 
Nycticebus. The new evidence derived from the pre- 
sent study is thus accordant with these evidences from 
previous studies described above. 

Ratajszczak (1998) proposed the theory to eluci- 


date the relationships between the two species, 
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